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ABSTRACT 

We consider a recently discovered class of instabilities, driven by cosmic ray streaming, in a variety of envi- 
ronments. We show that although these instabilities have been discussed primarily in the context of supernova 
driven interstellar shocks, they can also operate in the intergalactic medium and in galaxies with weak magnetic 
fields, where, as a strong source of helical magnetic fluctuations, they could contribute to the overall evolution 
of the magnetic field. Within the Milky Way, these instabilities are strongest in warm ionized gas, and appear 
to be weak in hot, low density gas unless the injection efficiency of cosmic rays is very high. 

Subject headings: cosmic rays — magnetic fields — instabilities 



1. INTRODUCTION 

Recently a powerful instability which couples a high flux of 
cosmic rays to their host medium has been discovered (Bell 
2004, Blasi & Amato 2008); we refer to this as the "Bell In- 
stability." The instability amplifies low frequency, right circu- 
larly polarized electromagnetic fluctuations with wavenumber 
parallel to the ambient magnetic field. In contrast to the clas- 
sical cyclotron resonant streaming instability (Wentzel 1968, 
Kulsrud & Pearce 1969), in which cosmic rays with Lorentz 
factor 7 amplify Alfven waves with wavelength of order the 
cosmic ray gyroradius r cr ~ 7c/cj„- (where cj„- is the non- 
relativistic ion-cyclotron frequency), the characteristic wave- 
length of the Bell instability is much less than r cr . The Bell 
instability is thought to be an important ingredient of diffu- 
sive shock acceleration in supernova remnants. It might be 
responsible for amplifying the magnetic field by up to ~ 2 or- 
ders of magnitude above its interstellar value and increasing 
the maximum energy to which cosmic rays can be accelerated, 
possibly up to the "knee" at ~ 10 15 eV (Drury 2005, Reville 
et al. 2008a). This is an important result, because it has been 
known since the work of Lagage & Cesarsky (1983) that stan- 
dard models of diffusive shock acceleration in the interstellar 
medium fail to reach the energy of the knee. 

Although the Bell instability has been applied primarily to 
supernova remnants, there are many other environments in 
which the cosmic ray flux may be large enough to excite it: 
galactic wind termination shocks, intergalactic shocks, and 
shocks in disks and jets. Cosmic rays streaming away from 
local sources or from their host galaxies may also constitute 
a sufficient flux. If the Bell instability exists in any of these 
environments it can amplify the magnetic field and transfer 
cosmic ray energy and momentum to the background plasma 
as well as increasing the efficiency of cosmic ray accelera- 
tion. In view of our current uncertainty as to how galactic and 
intergalactic magnetic fields originated and are maintained 
(Widrow 2002, Kulsrud & Zweibel 2008), any mechanism for 
amplifying them is of interest. Because fluctuations generated 
by the Bell instability have a definite helicity relative to the 
background magnetic field, they could be significant in am- 
plifying the field at large scales (Pouquet et al. 1976). This 
paper assesses the conditions under which the Bell instability, 
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or more generally any rapidly growing, nonresonant electro- 
magnetic streaming instability, can exist. 

In order to excite the Bell instability, the cosmic ray particle 
flux n cr vo, thermal ion density n,-, an d Alfv en speed va must 
satisfy the inequality (as we show in §12.2. U 

n cr v D >«,—-, (1) 
(7)c 

where (7) is of order the mean cosmic ray Lorentz factor. This 
can be written in terms of the cosmic ray and magnetic energy 
densities U cr , Ug, 



Equations ([T} and (01 express the requirement that the charac- 
teristic wavenumber of the Bell instability be much greater 
than the reciprocal of the mean cosmic ray gyroradius 

c(7)Mri. 

The Bell instability was originally derived for a cold 
plasma, which is valid when the thermal ion and electron gy- 
roradii r,, r e , are much less than the characteristic wavelength 
kseii of the instability; kseiir^e — > 0. Reville et al. (2008b) 
considered kr, small but nonzero. Their condition that the in- 
stability is significantly modified by thermal effects can be 
written in terms of the cosmic ray flux and ion thermal veloc- 
ity Vj = yflksT I nii are (as we derive in §I2.2.21 > 

v \ 

n cr v D >ni^r. (3) 
vf 

When the inequality (0) is satisfied, the wavenumber of max- 
imum instability, k W i ce (warm ions, cold electrons) decreases 
relative to kseii- And, while the maximum growth rate of the 
Bell instability, and the growth rate of the resonant streaming 
instability are independent of B as long as va much exceeds 
the Alfven speed Va, the growth rate of the thermally modified 
instability increases linearly with B. 

In this paper we extend Reville et al.'s analysis to cases 
where kr-, is not small, include cyclotron damping, and inves- 
tigate the properties of the instability by solving the plasma 
dispersion relation. The results are given schematically in 
Figure Q] and pr ecise ly for two representative environments, 
for (7) = 1, in §13.11 Equations ([J} and © define curves on 
the (n cr VD,B) plane. We show that these two curves, together 
with the requirements k wice r cr > 1, fc„, !CC r,- < 1, and k B eiiri < 1, 
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divide the (n cr VD,B) plane into the domains delineated in the 
figure. First, we consider the case of k W j ce r cr > 1: if the cos- 
mic ray flux is too low, resonant streaming instabilities can be 
excited, but nonresonant instabilities of the Bell or thermally 
modified Bell type are precluded. The threshold flux depends 
only on v, and the ion density n,, and is independent of B. 
Above this threshold, but below a second threshold, defined 
by k W i ce ri < 1, the plane is divided into three regions. If B 
is large enough that eqn. (fTJ is violated, nonresonant insta- 
bility is again precluded. Intermediate values of B, such that 
eqn. (01 is violated but eqn. (HJ is not, define the range of the 
standard Bell instability. If B is small enough that eqn. (01 
holds, there is nonresonant instability of the thermally modi- 
fied type. The condition k W j ce rj < 1 gives an upper limit to the 
flux at which nonresonant instability can be excited. Above 
this flux, the ions are unmagnetized and any instability which 
exists is driven solely by the electrons. This limit too, which 
is discussed further in §2 and in the Appendix, is independent 
of B and depends only on «, and v,. Evaluating these limits 
numerically, we find that there is a broad range of cosmic ray 
flux and magnetic fieldstrength for which the instability op- 
erates efficiently, should mediate the transfer of cosmic ray 
momentum and energy to the background, and should be a 
powerful source of circularly polarized magnetic field fluctu- 
ations. 

In §|2] we write down the general dispersion relation and 
solve it in various limits. In §[3] we give quantitative versions 
of FigureQ]and apply the results to several different astrophys- 
ical environments. Section [4] is a summary and discussion. 
Gaussian cgs units are used throughout. 

This paper is restricted to linear theory. The quasilinear 
and nonlinear evolution of the instability has been studied by 
a number of authors using both magnetohydrodynamic (Bell 
2004, 2005, Pelletier et al. 2006, Zirakashvili et al. 2008) and 
kinetic models (Niemiec et al. 2008, Riquelme & Spitkovsky 
2009, Ohira et al. 2009, Luo & Melrose 2009, Vladimirov et 
al. 2009). A host of saturation mechanisms have been investi- 
gated, including increase of the fluctuation scale to r~j, accel- 
eration of the thermal plasma, coupling to stable or damped 
modes, and relaxation of the exciting beam. Although un- 
derstanding saturation is essential for predicting the outcomes 
of the instability, delineating the regimes in parameter space 
where the instability exists is the first step. 

2. DERIVATION OF THE INSTABILITY 

We analyze the situation considered by Bell; a singly ion- 
ized plasma with ion number density «, and temperature T. 
There is a uniform magnetic field B = zB and a population of 
proton cosmic rays with number density n cr streaming along 
B with speed vp relative to the thermal ions. We assume from 
now on that the mean Lorentz factor (7) of the cosmic rays 
is order unity. There are applications, such as a "layered" 
shock precursor in which the most energetic particles have 
penetrated furthest upstream, where locally (7) ^> 1, and our 
results can be scaled readily to this situation; that application 
is important for determining the maximum energy to which 
particles can be accelerated. However, the instability will be 
excited by the bulk cosmic -ray population closer to the shock, 
which is important for field amplification. 

There are two populations of electrons, one with density «, 
which has no bulk velocity in the frame of the protons, and 
the other with density n cr which drifts with the cosmic rays at 
speed vo. Thus, the system is charge neutral and current free. 
This is the model used in Zweibel (2003) and Bell (2004). 



Amato & Blasi (2009) have considered the Bell instability 
when all the electrons drift at speed (n fr /n,)vo and found re- 
sults similar to those obtained for the two electron populations 
assumed here. But although the Bell instability in its original 
form is insensitive to the precise form of the thermal electron 
distribution function f e (v), it does turn out to depend on f e in 
a hot plasma. 

We briefly consider the constraints on f e in the Appendix. 
Based on an assessment of the electrostatic Langmuir insta- 
bility, we argue that for very large drifts and high cosmic 
ray densities, such as are expected in young supernova rem- 
nants expanding into diffuse interstellar gas, a separate elec- 
tron beam with density n cr and drift velocity Vd is unstable, 
and would tend to relax. For more moderate shocks, or in 
outflows where the drift speed is less than the electron thermal 
velocity, such a beam is stable. In general, stability consider- 
ations alone do not determine f e . 

The fastest growing nonresonant cosmic rays streaming in- 
stabilities are sensitive to the precise form of f e only in a 
plasma so hot and/or weakly magnetized that the thermal ions 
do not respond. Any instability present is then an instability of 
the electrons alone. Its physical significance is unclear, since 
it depends on the form of f e . 

2.1. Full Dispersion Relation 

We are interested in right circularly polarized electromag- 
netic fluctuations which propagate parallel to B and depend 
on z and t as expi(kz-uJt); thus instability corresponds to 
Im(cu) = LOj > 0. The dispersion relation for the fluctuations 
can be written as 
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where we have dropped the displacement current, oj pe j = 
(47r«,e 2 /m e .i) l l 2 are the thermal electron and ion plasma fre- 
quencies, us ce j = ^feB/m e i c are the electron and ion cyclotron 
frequencies, Z is the plasma dispersion function (Fried & 
Conte 1961) 
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(5) 



and the quantity £ r is defined in eqn. (A10) of (Zweibel 2003) 
and is plotted in Figure (|2). Although the exact behavior of 

depends on the cosmic ray distribution function, it has the 
general property that when kr cr 3> 1, C ~~ * ~1> with a small 
imaginary part of order (kr cr )~ l . The physics underlying this 
behavior is that cosmic rays barely respond to disturbances 
with wavelengths much less than their gyroradius, but the 
electrons do respond, resulting in a large perturbed current. 
A fraction (kr cr )~ l of cosmic rays have large enough pitch an- 
gles that they can resonate, resulting in a small imaginary part. 
As to the background plasma terms, the first term on the right 
hand side of eqn. (01 represents the response of the thermal 
ions, and the second term represents the cold electrons. 

We have solved eqn. (01 for a variety of ambient medium 
parameters and cosmic ray distribution functions, and have 
reproduced the plots of growth rate vs wavenumber in Bell 
(2004) and Blasi & Amato (2008). An example is shown in 
Figure0 which plots growth rate vs wavenumber at fixed cos- 
mic ray flux, ion density, and magnetic fieldstrength for three 
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FIG. 1. — A schematic of the instability regimes in the cosmic-ray flux vs. magnetic-field strength plane. At high magnetic-field strengths, in the unshaded 
region in this figure, the magnetic field is too strong to allow the non-resonant Bell Instability to grow. At lower magnetic-field strengths, a cosmic-ray flux 
dependent threshold (shown here with the solid line) is reached where the cosmic rays can launch the Bell instability; this limiting magnetic field scales as 
(Mcz-vd) 1 / 2 . At lower magnetic-field strengths (below the dashed line in this figure), thermal-pressure effects become important, and modify the maximum 
growth-rate of the Bell Instability; this threshold scales as (« c ,'Vo)' ■ All of these regions are bounded on the left and the right by limits on the cosmic-ray flux: 
to the left (kwicefcr > 1) is the limit that the cosmic-ray flux must be high enough to excite the non-resonant Bell Instability. To the right (i WICf r; < 1) is the limit 
beyond which the thermal ions are no longer magnetized; i.e., the gyroradius of the thermal ions is larger than the wavelength of the instability. The quantitative 
labels on this plane change with (7), T, and «,; quantitative results are given for particular physical cases, with {7) = 1, in Figures |6]and|2] 
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FIG. 2. — Imaginary part (top panel) and real part (bottom panel) of the 
function Q- for right circularly polarized waves, plotted vs dimensionless 
wavenumber kc/bJd, for the normalized cosmic ray distribution function 

<HP) = (4/"7>oXl +P 2 /PqT 2 , where p = m p c. 
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FIG. 3. — Growth rates Wj vs wavenumber k in units of r c j = u) c j/c 
for three different temperatures in a medium with B = 3fiG, m = lcm -3 , 
n„v D = 10 4 cm- 2 s" 1 . Solid line: T = 10 4 K, long dashed line: T = 10 6 K, 
short dashed line: T = 10 7 K. 



different temperatures: T = 10 4 K, T = 10 6 K, T = 10 7 K. The 
other parameters, B = 3 fxG, n, = 1 cm -3 , n cr vo = 10 4 cm -2 s -1 , 
are similar to those chosen by these authors as representative 
of cosmic ray acceleration in a young supernova remnant. 

The lowest temperature is essentially the cold plasma re- 
sult 3 . The growth rate peaks at kseii and then plunges, but 

3 At sufficiently low ISM temperatures ion-neutral damping of the fluctua- 
tions must also be considered; (Zweibel & Shull 1982, Reville et al. 2008b), 
but that is beyond the scope of this paper. 



remains positive as k increases. The T = 10 6 curve is simi- 
lar to the T = 10 4 K curve, except that the peak growth rate 
is reduced and occurs at a slightly smaller k, and the positive 
tail is damped. The T = 10 7 K curve is markedly different. 
The wavelength of the fastest growing mode is longer, and the 
peak growth rate is lower, than at T = 10 6 K. The instability 
cuts off abruptly at a longer wavelength than in the T = 10 4 K 
and T = 10 6 K cases, then re-emerges in a short interval of k 
before disappearing again. We explain these features with an 
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analytical treatment in the following two subsections. 

The effect of magnetic fieldstrength on the fastest growing 
mode is shown in Figures (@^,b). Figure (|4^) plots the maxi- 
mum growthrate LOugm vs B at T = 10 4 K, with other param- 
eters as in Figure (|3jl. For B < 1 /iG, uJijgm increases linearly 
with B. At larger B, u>ij g m is independent of B. Different be- 
havior is seen for kf gm : the wavenumber of the fastest growing 
mode, plotted in Figure (HJ)), is independent of B at low field- 
strength (here, for B < 1 /iG) and decreases as B~ l at larger B. 
These features, too, are derived in the next two subsections. 



The fastest growing mode occurs at the Bell wavenumber k Be u 
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FIG. 4. — Left panel: maximum growthrate of the streaming instability as a 
function of magnetic fieldstrength. Right panel: wavenumber k of the fastest 
growing mode. Here, T = 10 4 K and the other parameters are as in Figure[5] 



2.2. Analytical Results 
2.2.1. Standard case 

We recover the essential form of the Bell instability from 
eqn. (0]l by setting £ r = — 1 and taking the zero temperature 
limit of the right hand side. Using Z(z) — * -1 /z; \z\ S> 1, eqn. 
(0 becomes 



2 1-2 
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We now assume lu <C cj C! -, approximate (w ci +o;) 1 by w c /(l — 
ui/u> C i), and use u> 2 e /ui ce = -u) 2 pi juj C i. Multiplying the resulting 
dispersion relation by uo 2 v\/c 2 yields 
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The solution to eqn. ^} is 
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and has 



u(k Be u) = bJ BeU = 
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where the limiting expression on the last line holds for 

vd/va > 1. 

Equation ( TTOb shows that instability requires vo/va > 1. 
This is also the threshold for the classical resonant stream- 
ing instability. When the cosmic ray flux is too low to satisfy 
eqn. ([TJ, the dominant instabilities are resonant instabilities 
of Alfven waves. The peak growth rate occurs at wavenum- 
bers which resonate with cosmic rays near the mean energy, 
k ~ r~j , and is of order u> C i(n cr / «,)( Vd/ va - 1 ) (Kulsrud & Ce- 
sarsky 1971). For vd/va^> 1, this expression agrees to within 
a factor of order unity with u>Bell> but the two cannot be used 
simultaneously because they apply for opposite cases of the 
inequality (03 . 

At fluxes which satisfy eqn. (Q~|i, the assumptions made in 
deriving the classical resonant growth rate - that the under- 
lying waves are Alfven waves and that the growth time is 
much longer than the wave period - are incorrect (Zweibel 
1979, 2003, Achterberg 1983). At wavenumbers k < r" 1 , the 
growthrate is of order uj C i(n cr /ni)(vo/vA- 1) 1/<2 for vo slightly 
greater than v A and peaks at uj ci {n cr v D /nic) x l 2 for v D ^s> v A 
(Zweibel 2003). Comparing this expression to uiBeii, we see 
that u> res /ujBeU ~ { n i^Aj 1 n cr cvi,) x l 2 . From eqn. (Q]), we see that 
whenever the Bell instability operates, its growth rate exceeds 
the growth rate of the resonant instability, and the growth rate 
of the resonant instability is lower than predicted by the clas- 
sical theory. 

At k > 2kBeii, the nonresonant instability is stabilized by 
magnetic tension. Comparison with Figure [3] shows that al- 
though the growthrate plunges to low levels at this value of k, 
the instability does not completely disappear, as is predicted 
by eqn. (0. The oc k~ l tail seen in Figure[3]can be recov- 
ered by using the full ( r ; for kr cr ^> 1, £ r ~ — 1 + i/kr cr . 

The Bell instability is important in shock acceleration if 
the growth time is shorter than the time it takes the shock 
to travel through the layer within which the cosmic rays are 
confined. For acceleration at the maximum rate, the cosmic 
ray diffuse according to the Bohm formula with diffusion co- 
efficient D ~ cr cr , which sets the convection time across the 
layer as (uj C jvj) /c 2 )" 1 . The Bell instability will be able to grow 
if (n„ v D /«,v A ) > v 2 D /c 2 . 



2.2.2. Warm ions, cold electrons 

We now imagine decreasing B or increasing T such that 
kvi/\uj C i + bj\ is less than unity but not infinitesimal. In this 

case, Z(zd ~ — 1/z— l/2z 3 + i^Jne~~ . The imaginary part rep- 
resents ions for which the Doppler shifted wave frequency 
uj-kv = -ut C i\ i.e. ions in cyclotron resonance. The -l/2z 3 
term represents the finite gyroradius of the ions, which par- 
tially decouples them from the field. Using this approxima- 
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tion, eqn. © becomes 

2 / k 2 vf n cr . _o; 2 
V 2o7 C! - n,- £v; 

(11) 

Comparison of eqns. (Q and ( fTTT i shows that the thermal ions 
have two effects: cyclotron resonance, which is represented 
by the imaginary term, and a pressure-like effect due to the 
finite ion gyroradius, which increases the wave speed and is 
known as ion gyroviscosity. In the absence of cosmic rays, 
the dispersion relation agrees with Foote & Kulsrud (1979), 
and in the limit kr cr 3> 1, with Reville et al. (2008b), 4 except 
that ion cyclotron damping is neglected in both papers. 

Equation (fTTb represents the behavior of the growth rate 
quite well. This is shown in Figure (|5), which compares 07,- 
calculated from eqns. (0]i and ( fTTT i. 
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FIG. 5. — Left panel: Growth rate oj; vs scaled wavenumber kr cr as com- 
puted from the full dispersion relation eqn. J4} - long dashed line -, and the 
approximation jilt - short dashed line for T = 10 3 K and the other parameters 
as in eqn. {3). Right panel: same for T = 10 7 K. 



In the left panel, T = 10 3 K, thermal effects are unimpor- 
tant and the differences between the two curves are entirely 
due to the value adopted for C, r ; taking C, r = — 1 omits the ef- 
fect of resonant particles at high k and underestimates the cos- 
mic ray response at low k, thereby overestimating o;,. In the 
right panel, T = 10 7 K, thermal effects are important enough 
to modify the instability (see Figure [3] and below) but the an- 
alytical dispersion relation successfully reproduces the main 
features. Again, setting = -1 overestimates the growthrate 
at low k. 

It can be shown from eqn. (fTTT i that when 



< 



n cr v D 



1/3 



(12) 



V; V riiVi 

4 The sign of the second term on the left hand side of eqn. (3) in Reville et 
al. (2008b) is incorrect. 



the fastest growing mode is determined by competition be- 
tween the drift term 07„A;vo« c , /«, and the finite gyroradius 
term rather than the drift term and the magnetic tension term. 
This condition is equivalent to eqn. (01. The wavenumber 
k W ice of the fastest growing mode in this regime is 



k ■ 



n cr VP 
It 



1/3 



(13) 



Equation dT3l neglects cyclotron damping, and thus is valid 
only for k wice n - (n cr v D jnivi) 1 ^ < 1. As k wice n -> 1 cy- 
clotron damping becomes large, and the instability shuts off. 
A further constraint is that the instability be nonresonant: 
kwice^cr ~ (c/vi)k W i ce rj > 1. Both requirements can be writ- 
ten as temperatur e and density dependent limits on n cr vu', we 
return to them in 33.1l (eqns. [T7ll2"0"i i. 

The growth rate u W i ce corresponding to k wice is 



n cr vp\ 
m Vj J 



2/3 



(14) 



in agreement with Reville et al. (2008b). 

It can be shown from eqns. ( flOb . dT2l . and (fl4l ) that 
ijJwice/^Beii < 1- At the limits of validity of the warm ion ap- 
proximation, which is k W i ce Vi/u> c i = 1, uJ W i C e/<^Beli = Va/Vi- The 
suppression of the growth rate is due not to thermal ion cy- 
clotron damping, which is weak for bj ci /kv; 1, but to the 
restoring force exerted by the warm ions. Cyclotron damping 
does, however, come into play at shorter wavelength, obliter- 
ating the resonant tail of the instability. This happens roughly 

where (kvi/u C i)e~^"^ kv ^ > n cr vo/rijC. For the cosmic ray 
flux and ion density assumed in Figures [3] - |5j this occurs at 
kv-JiOci ~ 0.27, or kc/uj cj ~9x IQ 5 T~^ 2 . This is consistent 
with the behavior shown in Figure [3] 

At T = 10 7 K, something more complicated is going on: 
the instability growth rate decreases sharply after peaking 
near kc/u>d ~ 10 2 , as predicted by eqn. ( fT3l . but then has 
a brief resurgence. This is because the resonant ion cy- 
clotron term overwhelms the stabilizing gyroviscous term for 
Uci/kvi ~ 3/2, removing, in a small band of k space, gyrovis- 
cous stabilization. 

As we discussed in §2.2.1, the instability can only effi- 
ciently amplify magnetic fields at a shock if its growth time is 
faster than the convection time (u) c iV D /c 1 y~ . From eqn. ( TBI 
we see that this requires 



m Vi 



2/3 



(15) 



The criterion ( fT51 ) is independent of magnetic fieldstrength, 
and depends only on the cosmic ray flux, drift speed, and the 
density and temperature of the ambient medium. 

2.2.3. Hot ions 

When kr, > 1, eqn. ( fTTT i becomes invalid. In this limit, the 
argument of the plasma dispersion function in eqn. (0]l be- 
comes large, and Z(z) w -z~ l +i\f^. Physically, this means 
the ions are responding very little to the perturbation. Ion 
cyclotron damping is also weak, because the slope of the dis- 
tribution function is small at the resonant velocities. 

Under these conditions, the instability, if it exists at all, is 
due entirely to properties of the electron distribution function. 
In deriving the cosmic ray response function ( r used in eqn. 
(HI and plotted in Figure (O, we assumed the electrons are 
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cold and a fraction n cr /nj of them are drifting with the cosmic 
rays. As long as Vd/v c is sufficiently large, the electrons are 
unstable not only to the electromagnetic streaming instability 
considered here, but also to the much faster growing electro- 
static instabilities discussed in §|2] On the other hand, if the 
electrons were all drifting at speed n cr vo / V; both the electro- 
magnetic and electrostatic instabilities would be stabilized. 

Assuming the two peaked electron distribution, it can be 
shown that the instability growth rate is bounded above by 
u) C i(n cr Vo/niVi). This is generally less than uj w \ ce defined in 
eqn. (fT4l >. showing that these very short wavelength instabil- 
ities are not as important as the thermally modified or stan- 
dard Bell instabilities. At even shorter wavelengths, such that 
kr e > 1, the derivation of ( r becomes invalid. In view of our 
uncertainty about the electron distribution function, we have 
pursued the hot ion case no further. 

3. APPLICATIONS 
3.1. Instability regimes 

We begin by summarizing the different regimes of the 
streaming instability as functions of magnetic fieldstrength B 
and cosmic ray flux n cr Vu- These regimes were introduced 
without proof in |T]and depicted schematically in Figure (|T). 
The criteria used to delineate these regimes are approximate, 
but as Figures Hk,b indicate, the transitions between regimes 
are fairly sharp. 

According to eqn. dTJ, the condition that the Bell instability 
be nonresonant, i.e. that kB e iir cr > 1, is 

B<B s = 8.7x l(r 7 (n„.v D ) 1/2 , (16) 

where here and below n cr \'T> is given in units of cm" 2 s" 1 and B 
is in G. The condition that thermal effects modify the Bell in- 
stability such that the wavelength of the fastest growing mode 
is at k ~ ky/tce (eqn. |9]l rather than k ~ kseii (eqn. Qj) is 

B < Bm = 2.3 x 10"V /3 «! /6 («„.v D ) I/3 . (17) 

The condition for the thermally modified Bell instability to be 
nonresonant is k W j ce r cr > 1. At the same time, the thermal ions 
must be magnetized at k = k W i ce ; k W i ce ri < 1. These conditions 
limit n cr v D to the range 

v 4 

n~ <n cr v D <riiVi, (18) 
<r 

or numerically 

3 x KT'V-r 2 < n cr v D < 10 4 n,T 1/2 . (19) 

Finally, the condition that the ions be magnetized at k = kseii, 
kBeiin < 1, is n cr v D < HiV^/vi, or 

B>B T = 5x KT 10 r 1/4 (« c ,v D ) 1/2 . (20) 

Equations ( fT6b - ( f20b are plotted on the (n cr \>D, B) plane) in 
Figure s|6] and [7] Because eqns. ( fTTT i - d20l i depend on n t and 
T, we give two versions of the plot. Figure|6]represents the in- 
terstellar medium; n, = 1 cm" 3 , T = 10 4 K. Figure|7]represents 
the intracluster medium; n-, = 10~ 3 cm" 3 , T = 10 7 K. 

For n cr VD < n,-vf /c 3 , Bs < Bm- This is also the condition 
kwice r cr < 1 ( e q n - fl9t - and is represented by the leftmost ver- 
tical line on the plot. To the left of this line, there can be no 
nonresonant instability; since Bm > B$ the instability would 
take the thermally modified form, but for B > Bm the field is 



B 10" 6 
[Gauss] 10" 7 




10 -10 jq-7 1Q -4 1Q 2 1Q 5 1Q 8 

n cr vo [cm -2 s~ l ] 

FIG. 6. — Log-log plot of the (n cr vo,B) plane with vari ous regim es delin- 
eated. The quantities Bs, Bm, and Bt defined in eqns. j 161 11*71 . and 1201 
appear as lines with slope 1/2, 1/3, and 1/2, respectively. The flux limits de- 
fined in eqn. )18t appear as vertical lines. The plasma density is n, = 1 cm" 3 
and the temperature is T = 10 4 K. 
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FIG. 7. — Same as Figure|6]except «, = 10 3 cm 3 , T = 10 7 K 



too large and for B < Bm the flux is too low. Streaming insta- 
bility exists, but it is resonant. In the cold plasma limit, the 
Bell instability exists for any B < Bs, no matter how small. 

For n cr VD > rtivj/c 3 , Bm < Bs- For B > Bs, the field is too 
large for nonresonant instability (this is the case in the inter- 
stellar medium, away from cosmic ray sources). For B < Bs, 
the Bell instability operates in standard form as long as B ex- 
ceeds Bm and Bt- Although Bm < Bt is theoretically possible, 
it requires n cr v D > n^i, which is rather extreme. If we confine 
ourselves to n cr \'u < n(\>i (represented by the rightmost verti- 
cal line, defined by A; Hlce r, = lthen the Bell instability operates 
for Bm < B < Bs- For B < Bm, there is nonresonant instabil- 
ity as long as n cr va is to the left of the vertical line. To the 
right of this line, the ions are unmagnetized. As we have ar- 
gued, in this case the instability is controlled primarily by the 
electrons, and for fer, ~ 1, ion cyclotron damping is strong. 

In summary, nonresonant instabilities exist in the range of 
cosmic ray fluxes given by eqn. ( fl9] >. When B is between 
Bm, defined in eqn. dTvT >. and Bs, defined in eqn. ( TTBT ). the 
maximum growthrate is independent of B. When B < Bm the 
instability is thermally modified, occurs at longer wavelength, 
and grows at a rate proportional to B. This is also shown in 
Figure[8] which is a contour plot of the maximum growth rate 
in Equation (O as a function of cosmic -ray flux and magnetic 
field. Towards the lower-right of the plot, one can see the 
maximum growth rate decreasing, downward, with decreas- 
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ing B, but relatively constant above an approximately diag- 
onal line in (n cr vo,B) space from (10~ 10 cm" 2 s" 1 , 10~ 8 G) to 
(10" 4 ' 4 cm- 2 s-^KT^G). 




-8 -6 -4-2 2 
Logt/ic-Vfl/cnrV] 



FIG. 8. — Contour plot of the maximum growth rate (u)f gm ) of Eqn (4) 
plotted as a function of the logarithm of the cosmic-ray flux n c , i'o and the 
logarithm of the magnetic-field strength, B. The growth rate does not change 
with B in the region where the resonant and Bell Instability dominate, but 
at smaller magnetic-field strengths, u>f gm starts to decrease linearly with 
B. At relatively large magnetic field strengths and low cosmic-ray fluxes 
(B ~ l(r 6 G and n cr vo ~ 2.5 cm -2 s ), the growth rate falls slightly in the 
transition between the non-resonant Bell Instability and the resonant stream- 
ing instability, giving a slight curvature to the contours there. 



3.2. Astrophysical settings 

We now consider a few examples of specific astrophysical 
settings. In order to smoothly represent the transition from the 
standard Bell regime Bm < B < Bs to the thermally modified 
regime B < B M , we replace LU W j ce by Lo Be uB/B M , which agrees 
with eqn. (fTOt up to a factor of order unity. 

3.2.1. Supernova remnants and superbubbles 

First, we consider cosmic ray acceleration by a super- 
nova driven shock traveling at 10 4 km s" 1 through the in- 
terstellar medium. As in Figure [6^ we take n, = 1 cm -3 , 
T = 10 4 K. According to eqn. ( fl9] >, nonresonant instabil- 
ity exists for 3 x 10~ 8 cm~ 2 s _I < n cr \'r> < 10 6 cm" 2 s _I . Tak- 
ing n cr \'i) = 10 4 cm" 2 s _I , as assumed in Figure [3] we find 
from eqns. ( TTol l and dTvT > that the standard Bell instability 
operates for B between 1 . 1 and 87 /iG, which encompasses 
most of the fieldstrengths measured in the diffuse interstellar 
medium. The growth rate in the Bell regime is 2.5 x 10~ 4 s -1 
and 2.5 x 10~ 4 (fi/l.l ^G) in the thermally modified regime. 
The condition for efficient field amplification at shocks, u> > 
u> C i(vo/c) 2 , isuj> 1.1 x lO'B. Growth is amply fast in the Bell 
regime and also in the thermally modified Bell regime, since 
the growth rate and advection time scale in the same way with 
B. Of course, if B is too small, the cosmic ray acceleration 
time becomes long compared to the shock evolution time. But 
even if B ~ 10~ 9 G, about an order of magnitude less than the 



disordered field estimated by Rees (1987) to arise from a su- 
perposition of plerion supernova remnants in the early galaxy, 
the characteristic growth time of the instability is less than 
a year, much faster than the timescale on which the remnant 
evolves. 

On the other hand, if the shock propagates in a hot, low den- 
sity medium such as the superbubbles modelled by MacLow 
& McCray (1988), with n, = 3 x 10~ 3 cnT 3 , T = 3 x 10 6 K, 
then according to eqn. ( fl9] l, the nonresonant instability could 
exist for 8.1 x lO^cnrV 1 < n cr v D < 5.2 x K/cirrV 1 . If 
the cosmic ray injection efficiency were the same as at higher 
densities, n cr /n, = 10~ 5 , eqns. ( fTSI l and ( TT7l > show that the stan- 
dard Bell instability operates between 0.39 and 4.7 fiG. The 
upper limit is only slightly less than the rms Galactic field. 
The growth rate in the Bell regime, 1.4 x 10~ 5 , exceeds the 
advection rate through the cosmic ray scattering layer only for 
B < 1 3fiG, suggesting that nonresonant instability is less im- 
portant for shock acceleration in a low density medium than 
in a high density medium. The injection efficiency, rather than 
thermal effects, is the deciding factor: increasing n cr /tij above 
10~ 5 would enhance the growth rate in proportion. These re- 
sults imply that nonresonant instabilities might not occur for 
cosmic ray acceleration in superbubbles, where the combined 
effects of many hot star winds and explosions make the den- 
sity low. 

Denser superbubbles have been observed (Dunne et al. 
2001). Increasing «, to 3 x 10~ 2 cm" 3 while leaving the other 
parameters the same would increase lOBeii to 3.8 x 10~ 5 s" 1 
and change the range in which the standard Bell Instability 
operates to \.2\iG < B < 15/iG. The RMS Galactic field is 
estimated to be ~ 5.5/iG (Ferriere 2001), so the conditions 
for instability should be satisfied. The instability growth rate 
exceeds the advection rate through the acceleration layer as 
long as B < 4.1G. 

3.2.2. Shocks in galaxy clusters 

Next, we consider acceleration at shocks in the intraclus- 
ter medium. We take n, ~ 10~ 3 cm" 3 , T = 10 7 K. Accord- 
ing to eqn. ( fT9l , nonresonant instability exists for 3 x 
10~ 5 cm" 2 s _1 <rc cr V£><3.1 x 10 4 cm" 2 s _I . If acceleration oc- 
curs at the same efficiency assumed for galactic supernova 
remnants then n cr = 10~ 8 cm" 3 ; assuming vq ~ v$ ~ 1000 km 
s" 1 , then Bs, the maximum fieldstrength for nonresonant in- 
stability, is 8.7 x 10~ 7 G. This is slightly less than the fields 
inferred in galaxy clusters (e.g., Govoni & Feretti 2004), 
but the parameters are too uncertain to rule out resonant in- 
stability. According to eqn. ( flOl l, the growth rate in the 
Bell regime is 7.5 x 10~ 7 s" 1 . With these same parameters, 
Bm = 1 .6 x 10~ 7 G; below this value the instability is thermally 
modified and grows at the rate 7.5 x 10~ 7 (B/.16 /iG) s _1 . Due 
to the high temperature, the range in which the standard Bell 
instability operates without thermal effects is quite small, but 
the growth rate of the thermally modified instability is still 
fast. The condition for efficient field amplification at shocks, 
uj > L; ci (v D /c) 2 , is w > 1.1 X 10 _I B. Although the instability 
growth rates are lower than for supernova remnants, they are 
still large enough to satisfy this condition in both the standard 
and thermally modified nonresonant regimes. Thus, it appears 
that nonresonant instabilities could play a role in shock accel- 
eration, and could amplify magnetic fields, in galaxy clusters. 
This could make intergalactic shocks a favorable environment 
for acceleration of ultra high energy cosmic rays. 
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3.2.3. Unconfined galactic cosmic rays 

As our final example, we consider leakage of cosmic 
rays from galaxies into the intergalactic medium, which we 
take to have density and temperature n, = 10~ 6 cm" 3 , T = 
10 6 K (Richter et al. 2008). From eqn. CGD, nonreso- 
nant instabilities can be excited by cosmic ray fluxes be- 
tween 3 x 10~ 10 and 10 cm" 2 s" 1 . In the local interstellar 
medium, n cr ~ 10~ 9 cm" 3 , while vq is roughly the scale 
height of cosmic rays divided by their confinement time 
in the galaxy, or about 100 km s" 1 . This gives a galac- 
tic flux n cr vr> ~ 10~ 2 cm" 2 s _1 . Assuming cosmic rays em- 
anate isotropically from a characteristic galaxy size R g , we 
write n cr vo ~ 10~ 2 (R g /R) 2 (L cr /L crM w)- From eqns. dT6b 
and we find B s = 8.7 x lQ- s (R g /R)(L cr /L crMW ) 1 / 2 while 
B M = 7.3 x l0~ 9 (R g /R) 2/3 (L cr /L crMW ) 1 ' 3 . These values sug- 
gest that nonresonant instabilities could be excited in the in- 
tergalactic medium even if the fields are weaker than the 
10~ 9 - 10~ 10 G range often cited as upper limits (Kulsrud & 
Zweibel 2008). The growth rates, however, are rather slow. 
For example, if R/R g = 10, uJBeli = 2.5 x 10~ 9 \L cr /L cr Mw)- In 
this case, B M = 1.6 x 10" 9 G; for B = 10" 10 G, the maximum 
growth rate is about 4.1 x 10~ 10 '(L cr /L crM w) s_1 - Still, al- 
though the growth time exceeds 10 3 years, this is much shorter 
than any reasonable cosmic ray convection time. Therefore, 
nonresonant instabilities could be excited by cosmic rays from 
ordinary galaxies in the intergalactic medium at large. They 
could amplify intergalactic magnetic fields, and could heat the 
plasma. 

4. SUMMARY 

Nonresonant instability driven by cosmic ray streaming has 
emerged as a strong candidate for amplification of magnetic 
fields in environments such as strong shock waves, where the 
cosmic ray flux is large (Bell 2004). When the flux is high 
enough and/or the magnetic field is low enough, that eqn. 
(Q} is violated, the nonresonant instability replaces the clas- 
sical resonant streaming instability as the dominant electro- 
magnetic instability generated by cosmic rays. Although the 
instability scale length predicted by linear theory is small even 
compared to the cosmic ray gyroradius r cr , nonlinear simula- 
tions suggest that as the amplitude of the instability grows it 
generates fluctations at larger scales. This can increase the 
energy to which particles are accelerated in shocks. 

Cosmic ray acceleration and magnetic field growth are both 
of interest in a variety of environments, including young 
galaxies which may be actively forming stars but have not yet 
built up magnetic fields, shocks in galaxy clusters, and the in- 
tergalactic medium at large. In this paper we have carried out 
a parameter study of nonresonant instabilities including ion 
thermal effects. We solved the full dispersion relation © nu- 



merically and verified that a simple analytical approximation, 
eqn. (fTTl i. is quite accurate in the wavenumber regime of in- 
terest. We corroborated the criterion of Reville et al. (2008b) 
for when ion gyroviscosity reduces the instability growth rate 
and shifts it to longer wavelength. We showed that ion cy- 
clotron damping cuts off the instability at short wavelengths 
and argued that at wavelengths short enough that the ions are 
unmagnetized the instability depends only on the electron dis- 
tribution function, the prediction of which is beyond the scope 
of this paper. 

The joint requirements that the instability wavelength be 
much less than the cosmic ray gyroradius but much more than 
the thermal ion gyroradius limits the range of fluxes which 
excite nonresonant instability to n,vf/c 3 < n cr vo < In 
practical terms, this range is large and accommodates most 
cases of interest. Within the unstable range, there is a "strong 
field" regime in which all streaming instabilities are resonant, 
an "intermediate" regime in which nonresonant instability in 
the form derived by Bell dominates, and a "weak field" regime 
in which the instability is thermally modified. In the Bell 
regime the maximum growth rate is independent of B but 
in the thermally modified regime it depends linearly on B. 
Young galactic supernova remnants are generally in the inter- 
mediate regime unless the ambient medium is hot and rarefied 
(like the interior of a superbubble), in which case the instabil- 
ity is weakened. Generally, if B is in the nanogauss range, the 
growthrates are fast enough for nonresonant instability to be 
a potential source of magnetic field amplification in weakly 
magnetized interstellar and intergalactic gas. At much lower 
fieldstrengths, the instability is too slow to be of interest, but 
other instabilities, such as Weibel modes, could be an impor- 
tant ingredient in magnetogenesis (Medvedev et al. 2006). 

Although nonresonant instabilities amplify magnetic fields 
on rather small scales - much smaller than the eddy scales 
characteristic of interstellar and intergalactic turbulence - they 
should not be ignored in discussions of magnetogenesis. Be- 
cause only the right circularly polarized modes are unstable, 
nonresonant instabilities are a source of magnetic helicity on 
scales at which the background magnetic field is coherent. 
Magnetic helicity is thought to be a key ingredient in the 
growth of large scale magnetic fields from small scale fluc- 
tuations (Pouquet et al. 1976). Cosmic ray generated fluc- 
tuations could be important in driving an inverse cascade of 
magnetic power to longer wavelengths, and could prevent the 
pileup of power at short wavelengths that currently confounds 
interstellar and intergalactic dynamo theories. 
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APPENDIX 

APPENDIX: ELECTRON DISTRIBUTION FUNCTION 

Here we briefly consider the constraints on the electron distribution function f e 

One way or another, the cosmic ray current must be cancelled: an uncompensated cosmic ray current en cr vo flowing in a 
channel of width L pc measured in parsecs generates a magnetic field B ~ Q.5n cr v D L pc G. Even the galactic flux of 10~ 2 cm -2 s -1 
with L pc = 1 would generate a 5 mG field. Since cosmic rays are ion dominated, thermal electrons must cancel their flux. 

When an electron beam drifts with respect to the bulk plasma, it can excite rapidly growing electrostatic instabilities which 
tend to redistribute electron momentum and bring the system to a state of marginal stability. Langmuir waves (also called plasma 
oscillations) with wavenumber ui pe /v (u> pe is the electron plasma frequency (Ai:n e e 2 /m e ) 1 ' 2 ) are destabilized if df e /dv > 0. If 
the beam and bulk electrons have the same temperature T e and the beam velocity Vb much exceeds the electron thermal velocity 
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v e = w2kjjTe/m e (which is necessary for instability if the beam density rt/, is much less than the bulk density n e , the case of 
interest here), then the requirement for stability is approximately 



V ±e*l* > "± (Al) 



II,, 



(e.g. Krall & Trivelpiece 1973). In shock acceleration, it is sometimes assumed tib/n e ~ 10 5 , According to eqn. (lAU . stability 
then requires Vb/v e < 3.5. Assuming T e = 7) in the upstream plasma, beams associated with shocks of Mach number M < 

2>.5\Jrrii/m e are stable while shocks at higher M are unstable. In a 10 4 K gas, the stability boundary is at about 1500 km s -1 . Thus, 
while the Langmuir instability is a constraint for very fast shocks, it is probably irrelevant for older supernova remnants, and in 
galaxy cluster accretion shocks or galactic wind termination shocks, where the background gas is hot and the Mach numbers are 
expected to be moderate. The fluxes associated with cosmic ray escape from galaxies are probably also electrostatically stable. 

Therefore, it appears that f e is not determined by stability considerations alone, but depends on other factors such as the history 
of the system and the source of cosmic rays. 
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